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1. Pointing Accuracy of the Low-Cost Large Space Telescope Due to 
Noise and Quantization 

1 - 1 . Introduction 

The objective of this study is to conduct an investigation on the 
pointing stability of the low-cost Large Space Telescope (LST) system. 

The low-cost LST is characterized by the use of reaction wheels for 
the generation of control torques. Because of the critical requirement 
on the pointing accuracy of the LST, the nonlinear frictional characteristics 
of the bearings of the reaction wheels cannot be neglected. It is well 
known that the nonlinear friction can cause limit cycles in a closed-loop 
system. 

Another possible source of pointing error in the LST is due to the 
effect of quantization and sensor noise. Since the LST is a digital 
system, D-A and A-D converters, and sensors for positional and rate 
feedbacks are used. Sensor noise and amplitude quantization will also 
cause pointing error in the LST. In addition, quantization is a 
nonlinear phenomenon so that it may also cause self-sustained oscillations 
in the closed-loop system. 

The dynamic modeling of the single-axis LST is described in this 
chapter. Several methods of evaluating the attitude error of the digital 
LST due to quantization and noise inputs are given. 
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1-2. Dynamic Models of the Low-Cost LST System 


The dynamic model of the single-axis low-cost LST system with 
sampled data is shown in the block diagram of Fig. 1-1. The rigid 
body is represented by the double-integrator transfer function. The 
controller is formed by proportional , rate, and integral feedbacks of 
the vehicle attitude. The nonlinear element N in the reaction wheel 
dynamics represents the rolling friction, and its functional description 
is given by the well-known Dahl model. 

The definitions of variables and the values of parameters and 
constants are tabulated) in Table 1-1 [!]. 


^c 

4> b 

4> b 

T 

c 

T RW 

0 RW 

0 RW 

T F 

^e 

“p 

K R 


Table 1-1. 

Reference input command 
Body attitude of LST 
Body rate of LST 
Torque command of reaction wheel 
Torque output of reaction wheel 
Angular velocity of reaction wheel 
Angular displacement of reaction wheel 
Frictional torque of reaction wheel 

1 .65 x 10^ N-m/rad 
3.71 x 10 5 N-m/rad/sec 
7.33 x 10^ N-m/rad/sec 
0.484 ohms 
10000 


Attitude error 

Proportional gain of controller 
Rate gain of controller 
Integral gain of controller 
Feedback resistance 
Voltage amplifier gain 
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4 


“M 

k b 

J RW 


Motor torque constant 
Back emf constant 

Moment of inertia of reaction wheel 

Moment of inertia of vehicle about 
pitch axis 

Motor time constant 

Motor resistance 


0.484 N-m/amp 
0.484 volt/rad/sec 
0.2 Kg-m 2 

41822 Kg-m 2 
0.002 sec 
10 ohms 


Although the low-cost LST system has digital control, it is 
informative to analyze the system of Fig. 1-1 first without the sample- 
and-hold device. Figure 1-2 shows the signal flow graph of the continuous- 
data LST system. The characteristic equation of the system is determined 
from Fig. 1-2, 


where 


Wk 

A - 1 + RpG 2 G 3 + K b K m G 3 G 6 + ^ G 2 G 3 G 4 + R F G 1 G 2 G 3 G 4 G 5 +N ^ G 6 G 7 + R F G 2 G 3 G 6 G 7^ “ 0 

(1-D 


G] = K p + Kj/s 
G 2 = K a /s 


> 3 R + (l + x s ) 
v m 


G 4 Js 
v 


V? 


G^ = 


6 J RW S 


G 7 ■? 


and N denotes the analog describing function of the reaction wheel nonlinearity. 
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For the linear model of the LST system, we set N = 0, and the 
characteristic equation becomes 


VW 5 + W 4 + < J v R F K A + + r fV/ 

+ k a r f k p s + k a k i r f ■ 0 


(1-2) 


It is of interest to investigate the dynamics' of the reaction wheel 
and the vehicle. The open-loop transfer function between T c and <j>g is 




R F K A 


V 2 ( r tV 2 + r t s + ¥a + ^ 


(1-3) 


J RW 


Substituting the system parameters into Eq. (1-3) gives 


= 5.7865 

T c s 2 (s 2 + 500s + 242058.5) 

5.7865 

s 2 (s + 250 + j423.74) (s + 250 + j423.74) 


(1-4) 


Thus, the reaction wheel is shown to have relatively fast dynamics. 

Substitution of the values of the system parameters into Eq. (1-2), 
and simplifying, the characteristic equation of the linear LST system' is 
written 


s 5 + 500s 4 + 242058.5s 3 + 2.1468 x 10 6 s 2 + 9.5476 x 10 6 s 

+ 4.24145 x 10 6 = 0 


0-5) 


The roots of the characteristic equation are: 
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s = -0.49659 
s = -4.22743 + j‘4.25123 
s = -4.22743 - j4. 25123 
s = -245.524 + j421. 118 
s = -245.524 - J421.118 

Note that the damping ratio of the dominant complex roots is 
0.705, and the natural undamped frequency is 6 rad/sec or 0.954 Hz. 

These parameters are achieved by selecting the controller constants K p , 

K r and Kj at the indicated values. However, the poles of the reaction 
wheel dynamics at s = -250 + j423.74 and s = -250 - j423.74 are only 
slightly affected by the body controller and they account for the charac- - 
teristic roots at s = -245.524 + j 421 .118 and s = -245.524 - J421.118 of 
the overall closed- loop system. Since these fast roots are very far away 
from the dominant ones, this means that for all practical purposes the 
dynamics of the reaction wheel can be neglected as far as the linear 
system is concerned. Figure 1-3 shows the block diagram of the simplified 
continuous-data low-cost LST system, and the digital system is shown in 
Fig. 1-4. 

The closed-loop transfer function of the continuous-data system of 
Fig. 1-3 is 


V 3 * 



V + K i 

+ k r j2 + v + K I 


(1-6) 


For the digital system of Fig. 1-4, the closed-loop transfer function 


is written 
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$ B (z) 6 a (z)G b (z) 

$ c (z) “ 1 + G c (z) + G a (z)G b (z) 


(1-7) 


where 


G ho< 5 > ■ 


-Ts 


g a ( 2) ’i (6 ho G i> = (1 - z '’>3 


KjS + Kj 


K z + K t T - K 

_e I e. 

z - 1 


(1-8) 


G B^ ^ G ho G 4 G 5^ ( ] " z )J 


f >. 'k 


1 


J s 3 
l V J 


= r(z + D 

2J v (z - 1) : 


(1-9) 


G c< z > =$< G ho G 4 K R> ■ K rC - Z '’>J 


r ^ 


V 


"R T 


J v z - 1 


(1-10) 


Equation (1-7) is simplified to 


V Z) 


T 2 (K p z 2 + KjTz + KjT - K p ) 


<1> c (z) 2J v z 3 + (T 2 K p + 2K r T - 6J v )z 2 + (6J y - 4K R T + T 2 K I T)z + (2K R T + KjT 3 - 2J y - K p T 2 ) 


( 1 - 11 ) 


The characteristic equation of the system is 


2J v z 3 + (T 2 K p + 2K r T - 6J v )z 2 + (6J y - 4K R T + T 3 Kj)z 


+ 2K r T + KjT 3 - 2J v - K p T 2 =0 (1-12) 


Or, 
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83644z 3 + (1.65 x 10 6 T 2 + 7.42 x 1 0 5 T - 2.50932 x 10 5 )z 2 

+ (2.50932 x 10 5 - 14.84 x 10 5 T + 7.33 x 10 5 T 3 )z 

+ (-8.3644 x 10 4 + 7.42 x ]0 5 T + 7.33 x 70 5 T 3 - 1.65 x 10 6 T 2 ) = 0 

(1-13) 

The characteristic equation roots are tabulated below as functions of 
the sampling period T: 


T (msec) 


Roots 

0.1 

1, 

0.998 ± jO.0025 

0.5 

1, 

0.998 ± jO. 0022 

1 

1, 

0.996 ± jO.0043 

5 

0.998, 

0.979 ± jO.0208 

10 

0.995, 

0.957 ± j 0.041 1 

25 

0.987, 

0.889 ± jO. 0988 

50 

0.975, 

0.766 ± jO. 1825 

100 

0.950, 

0.483 ± jO. 2901 

120 

0.940, 

0.356 ± j 0.300 

150 

0.925, 

0.150 ± jO. 260 

170 

0.915, 

0.0034 ± jO. 144 

180 

0.910, 

0.049, -0.195 

200 

0.900, 

0.107, -0.571 

220 

0.890, 

0.1116, -0.9078 

230 

0.885, 

0.1085, -1.077 (unstable) 


The root loci of Eq. (1-13) are sketched in Fig. 1-5 as a function of 
T. For small sampling periods, the characteristic roots are all located near 

















the z = 1 point. The linear digital LST system becomes unstable when 
T exceeds approximately 225 msec. 



13 


1-3. Effects of Quantization on Pointing Stability of the Low-Cost LST - 
Limit Cycle Conditions 

Quantization occurs at at least three places in the LST system. Two 
are at the displacement sensor and the rate sensor where A-to-D converters 
are used. A quantizer is also needed at the control torque input to the 
reaction wheel since D-to-A conversion is effected there. In addition, if 
the integral control Kj/s is implemented digitally, quantization should be 
considered in the digital controller as well. Figure 1-6 shows the digital 
LST system with quantizers. The z-transfer functions are defined in 
Eqs. (1-8), (1-9), and (1-10). The quantizers at the displacement 
sensor, the rate sensor, and the reaction wheel control torque are denoted 
by Qp, Q R , and Q-p respectively. The input-output relation of a quantizer 
is shown in Fig. 1-7. The quantization level is represented by h. 

We shall analyze the effects of quantization on pointing stability 
or accuracy of the LST by means of three different methods. The first 
method utilizes the deterministic approach and establishes a least upper 
bound on the pointing error due to quantization. The second method relies 
on treating the quantizer as a noise source, and statistical analysis is 
applied. The third method is also a statistical approach which represents 
the quantizer by a linearized gain K e ^(z). 

It should be pointed out that a system with quantization is a nonlinear 
system and its behavior cannot be predicted by linear theory. One of the 
well-known phenomena of a nonlinear system is that sustained oscillations 
may occur. When a digital system has several quantizers, it is extremely 
difficult to predict the condition of self-sustained oscillations. 
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To illustrate the effects of quantization, and how quantization 
can cause sustained oscillations in an otherwise stable linear system, 
let us refer to the digital systems shown in Figs. 1-8 and 1-9. The 
difference between the two systems in Figs. l-8a and l-9a is that the 
former has negative feedback and the latter has positive feedback; 
but both systems are stable. 

For r(k) = 0, both linear systems have zero steady-state values 
for c(k); that is, c(k) = 0 for k-*», for arbitrary initial state 
c(0). He shall show that when quantization is considered, the system 
in Fig. 1 -9b has a steady-state error, whereas the system in Fig. 1 -8b 
exhibits a sustained oscillation. 

Let the quantization level h be 2 in Fig. 1-7, and when the input 
is an integer, the output is the same integer. The state equation of 
the system in Fig. 1 -9b is 

c(k+l ) = Q[0. 9c(k)] (1-14) 

For c(0) = 10, it can be easily shown that c(k) = 4 for k _> 2. 

The state equation of the system in Fig. 1 -8b is 

c(k+l ) = Q[-0.9c(k)] (1-15) 

For c(0) = 10, c(k) = 8 for k _> 2 even, and c(k) = -8 for k > 2 
odd. Thus, the state of the system oscillates between -8 and +8. 
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1-4. Effects of Quantization on Pointing Stability of the Digital Low- 
Cost LST - Least Upper Bound Quantization Error 

In this section the effects of quantization on the low-cost LST 
system are investigated using a deterministic approach. The method 
of analysis is based on the "worst" error condition due to quantization. 

In general, the analysis gives a conservative estimation of the quantization 
error. 

Since the quantization error has a maximum bound of ±h/2, the "worst" 
error due to quantization in a digital system can be studied by replacing 
the quantizer in the state diagram by a branch with unity gain and an 
external noise source with a signal magnitude of ±h/2. The block diagram 
of the digital LST system with quantizers shown in Fig. 1-6 is redrawn 
in Fig. 1-10 with the noise sources. The transfer functions G A (z), Gg(z), 
and G c (z) are defined in Eqs. (1-8), (1-9), and (1-10), respectively. 

The z-transform of the body attitude of the LST due to the three 
quantizers when <J> c = 0 is 


hf, h- hr, 

z 4 T G A^ 1 T 1 ~2 
V z) = 1 + G a (z)G b (z) + G c (z)' G B ( z) 


Substitution of Eqs. (1-8) through (1-10) into Eq. (1-16), and 
simplifying, we have 


$ B (z) 


z [±^(K d z + K i T- K p )i (^+^)(z-1)]T 2 (z+1) 
Z " 1 2J v (z-l) 3 + T 2 (z+l)(K p z + K I T- K p )+2TK R (z-l) 2 


(1-16) 


(1-17) 
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The steady-state error of the body attitude due to the quantization effects 
is obtained by applying the final-value theorem to Eq. (1-17), (if the 
system is stable). 

Thus, 

- 1 

tun $ g (kT) = tim (1 - z )$g(z) = ±-p~ (1-18) 

k-*» z-*l 

It is interesting to explore the significance of this result on 
the error due to quantization. Firstly, the quantization error at the 
displacement sensor is propagated through the system without change in 
amplitude. Secondly, the errors due to the torque and rate sensor 
quantizers are completely eliminated at the output position. This is 
attributed to the integral control Kj/s in the forward path. 

Vtgttal Impl&myitcution ofi Fonmnd Contn.oM.eA 


If the proportional -plus-integral controller is implemented 
digitally, the transfer function G^(z) becomes 


G a (z) = K p 


TKt ( z + 1) 
+ 2(z - 1) 


(1-19) 


Substituting G A (z) from Eq. (1-19) into Eq. (1-16), we can again show 
that the steady-state error in $g(kT) due to the three quantizers is ' 
±hp/2, and the errors due to and Q T are completely eliminated. 

In reality, the digital implementation of the controller should 
also include quantization in the digital process. Figure 1-11 shows 
the block diagram of the LST system with quantizations also considered 
in the digital controller. The z-transform of the body attitude of the 
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system with quantizers in the digital controller 
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LST due to the five quantizers when $ c = 0 is 



( 1 - 20 ) 


where G^(z) is given by Eq. (1-19) and Gg(z) and G c (z) are given in 
Eqs. (1-9) and (1-10), respectively. 

Applying the final-value theorem to Eq. (1-20), we have 


£-on ( kT ) = 

k- B 


+ h 2 + hpKjT 



(1-21) 


It is interesting to note that the introduction of the integral 
control eliminates the noise signals that enter at all points after the 
integral control in the control loop; however, the digital implementation 
of this control in turn produces a quantization error which is bounded 
by ±h 2 /2KjT, where h 2 is the quantization level. 
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1-5. Effects of Quantization on Pointing Stability of the Digital Low- 
Cost LST - Equivalent Noise Source, Statistical Analysis 

In this section the pointing errors of the digital LST system due 
to quantizations in the displacement, rate, and torque channels are 
investigated by statistical means. The rms (root-mean-square) error 
in 4>g due to quantizations that are represented by equivalent Gaussian 
noise sources with zero mean, or white noise sources, is determined by 
setting <)> = 0. The results are then compared with those of the continuous- 

data LST obtained in reference [1]. 

The block diagram model of the continuous-data LST with the quantizers 
replaced by equivalent noise sources is shown in Fig. 1-12. The equivalent 
digital system is shown in Fig. 1-13. It is assumed that the equivalent 
noise sources that represent the quantization operations are white, so 
that their power spectral density functions in both the s and the z 
domains are constants. Therefore, 

$p(s), $q(z) = = power spectral density of displacement quantizer 

<^(s), <I> R (z) = = power spectral density of rate quantizer 

$y(s) , $y(z) = $y = power spectral density of rate quantizer 


Let the rms attitude error of the continuous-data LST due to 
be represented by Og^. The rms errors are given by the following 
relations: 

Displacement : 


BD 


2JJ f . M 0 (s)M D (-s)ds • * D 

-J°° 


1/2 


( 1 - 22 ) 
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where M^(s) denotes the transfer function between $q(s) and $ B (s), 


* R (S) 




V + K R = 2 + V + K I 


0-23). 


Similarly, 

Rate: 


BR 


2irj 


-J c 


M R ( S )M r (~ s )d s • $ F 


1/2 


0-24) 


Torque: 


f 


BT 


1 

2uj 


r J c 


'J°° 


Mj(s)M T (-s)ds • 


1/2 


(1-25) 


where 


■Mr< S ) 


- Ms] = 


<M S ) 


■T'^ , -3 . .2 


v + y + y + k i 


0-26) 


This result is obtained with the assumption that the quantizer noise 
is not multiplied by the sensor gain K R , as in reference [1]; otherwise, 
the right-hand side of Eq. (1-26) should be multiplied by K R . 

The total rms attitude error due to all three quantizers is simply 
the sum of the errors due to each noise source acting alone; that is. 


°B " °BD + a BR + g BT 


0-27) 


From Eq. (1-26) it is easy to see that a BR and a B j have the same magnitude. 
The line integral of the form of Eq. (1-22) can be evaluated by 
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contour integration, and then the residue theorem. An integral table is 
also available for the evaluation of the contour integral. Using Eqs. 
(1-23) and (1-26), the rms attitude errors of the continuous-data system 
are obtained as 


BO 


K 2 + 

^4“ 


K I K R 


, 1/2 


W 


(2.5748£ d ) 1/2 = 1.605 /$ d 


(1-28) 



1 




1/2 


= (0.8598 x 10-\) 1/2 
= 0.927 x 10~ 6 /^ 


(1-29) 


a BR = 0.927 x ]Q- b J¥ R 


(1-30) 


For the digital LST, the rms attitude errors are: 
Displacement : 


BD 


2M M D (z)M D (z- 1 )z _1 dz • 4 d 


1/2 


(1-31) 


Rate: 


T* - 

7 br 


[ Sj|V z)M R (2 ’ 1)z ‘ ldz ‘ *R 


1/2 


(1-32) 


Torque: 


T* = 

r BT 


Sjj M T (z)M T U-’)^ 1 dz • $ T 


1/2 


0-33) 


The transfer functions, M^(z), M R (z), and My(z) are determined from 
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Fig. 1-13. 


M d (z) 


G a (z)G b (z) 

1 + G a (z)G b (z) + G c (z) 


(1-34) 


g b (z ) 

-M r (z) = M t (z) = i + g a (z)G b (z) + G c (z) 


(1-35) 


where 


g b ( z) 


= T z (z + 1) 
2J v (z - I) 2 


(1-36) 


G c (z > 


K rT 


(1-37) 


and 


G fl (z) 


V ♦ K,T - K p 
z - 1 


(sample-and-hold and 
analog controller) 


(1-38) 


(2K + TKjz + TKj - 2K 

G fl (z) = y ? / _ ■ -py ^ (digital implementation (1-39) 

1 v - ; 0 f controller) 

Substituting Eqs. (1-36) through (1-39) into Eqs. (1-34) and (1-35), 
we get the following transfer functions which are used for the computation 
of the rms attitute errors in Eqs. (1-31) through (1-33). , , 

With zero-order hold and analog controller. 


T 2 [K z 2 + K t z + (K T - KJ] 
M ( z ) = 2 — 2 

V z; A 1 


(1-40) 
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¥ z > ■ V z) - J - ( 2 a ) ' - U 


0-41) 


where 


A 1 = 2J v z 3 + (-6J V + 2K r T + T 2 K p )z 2 + (6J v - 4K R T + KjT 3 )z 
+ (-2J v + 2K R T + KjT 3 - K p T 2 ) 


(1-42) 


With digitally implemented controller. 


V z) 


T 2 [(2K p + TKj )z 2 + 2TKjZ + TKj - 2K p 

V 


(1-43) 


•M r (z) = M t (z) = 2jZ - { \ — U 


(1-44) 


where 

*2 = 4J v z 3 + ( - 1 2 J v + 2T 2 K p + T 3 Kj + 4K r T)z 2 

+ (1 2 J + 2T 3 K t - 8K d T)z + (4K d T + T 3 K t - 2T 2 K - 4JJ 

V IK K 1 p V 

(1-45) 

A tabulation method or a numerical method [3] can be used to 
evaluate the contour integrals of Eqs. (1-31), (1-32), and (1-33), once 
the values of the parameters of the transfer functions are known. Table 

1-2 gives the values of ( <f> B ^ (^msH and <*Brms ) R for var1ous 
values of the sampling period T. 
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Table 1-2. 

T Sample-and-Hold with Digitally Implemented 

(msec) Analog Controller Controller 



°BD 

O* 

°BT 


°BD 

a* 

°BT 


0.01 

0.0052 

2.932 

X 

10" 9 

0.0052 

2.932 

X 

10“ 9 

0.1 

0.0161 

9.274 

X 

10‘ 9 

0.0161 

9.274 

X 

10' 9 

1.0 

0.0508 

2.936 

X 

10- 8 

0.0508 

2.936 

X 

10" 8 

5.0 

0.1140 

6.599 

X 

10- 8 

0.1141 

6.595 

X 

10 -8 

10 

0.1621 

9.394 

X 

10' 8 

0.1622 

9.383 

X 

10" 8 

25 

0.2602 

1.515 

X 

o 

1 

0.2607 

1.511 

X 

10“ 7 

50 

0.3777 

2.216 

X 

o 

1 

0.3796 

2.205 

X 

10 -7 

100 

0.5652 

3.368 

X 

10~ 7 

0.5720 

3.339 

X 

10" 7 


The results tabulated in Table 1-2 show that the rms errors of the 
LST with sample-and-hol d and the analog controller are very close to those 
of the LST with the controller implemented digitally. 

It is interesting to show that the rms attitude errors of the digital 
LST due to quantization are related to those of the continuous-data system 
in Fig . 1-12. 

Applying the limit as T approaches zero to Eq. (1-31), we get 


ZMn cr* n = tlm 
T-0 BD T+0 



M D (z)M D (z“ 1 )z' 1 dz$ D 


1/2 


Since 


Um M n (z) = M n (s) 
T-*0 U U 


(1-46) 


(1-47) 
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where Mp(s) is given in Eq. (1-23), and Mp(z) is given by either Eq. 
(1-40) or (1-43), depending on. the way the controller is implemented, 
and 


z ^dz = T ds 


(1-48) 


Eq. (1-46) is written 


ton a* 
T-*0 BD 


r J 00 


2rrj J- j 


M D (s)M Q (-s)ds$ 


1/2 


= A o 


BD 


(1-49) 


or 


'BD 


-= tim a* 

/r T _ 0 BD 


(1-50) 


The meaning of this relation is that the mean-square value of the 
attitude error of the continuous-data system is equal to 1/T times the 
attitude error of the digital system as T approaches zero. Table 1-3 


gives the values of a*p//T$D and for various values of T. 

These values are very close to the values of Ogp = 1.605 and Ogy = 

- fi 

0.926 x io , respectively, especially at very small sampling periods. 

These relations show that, in general, for small sampling periods, 
the attitude error due to a white noise input in a digital system will 
be less than that of the same system without sampling. For example, 
for T = 25 msec, the rms attitude error is 0.2602/$^ , whereas 
Ogp is 1.605/4^ . Therefore, the use of the continuous-data LST model 
of Fig. 1-12, as in [1], for the error analysis of the digital LST 
results in conservative results. 
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Table 1-3. 


T 

(msec) 

a BD 

At 0 

a* 

°BT 

/T4j- 

0.01 

1.644 

0.927 

x 10" 6 

0.1 

1.610 

0.927 

x 10" 6 

1.0 

1.606 

0.928 

x 10 -6 

5.0 

1 . 61 2 

0.933 

x 10" 6 

10 

1.621 

0.939 

x 10" 6 

25 

1.645 

0.958 

x 10' 6 

50 

1.689 

0.991 

x 10" 6 

100 

1.787 

1.065 

x 10 -6 


The results in Tables 1-2 and 1-3 again show that for the same 
quantization levels, the quantizer in the displacement channel produces 
far greater attitude errors than those due to the quantizers in the rate 
and torque channels. 
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1-6. Effects of Quantization on Pointing Stability of the Digital Low- 
Cost LST - Quasi! inear Analysis, Statistical 


In this section the effects of quantization on the low-cost LST 
are studied by means of a quasi 1 inearized equivalent-gain approach [4], 

[5], Although the analysis is conducted in the statistical sense, there 
is a basic difference between the present analysis and the one conducted 
in Section 1-5. In the previous section, the rms attitude error of the 
LST is evaluated by treating the quantizers as noise sources that are 
stationary Gaussian processes with zero means. In this section, 
the quantizer will be treated as a nonlinear element whose input is 
a random stochastic process. An equivalent gain, K (z), is derived for 
the quantizer. The attitude of the LST is then determined with the quantizer 
replaced by K gq (z) , and when the system is subject to a stochastic input. 

The only restriction with this method is that only one quantizer can be 
considered at a time. Since it has been established that the attitude 
of the LST is more sensitive to the quantizer Q^, we shall consider only 
the system model shown in Fig. 1-14. 

The quantizer Q q is isolated as shown by the block diagram in 
Fig. 1 - 1 5a ; the input is x(t) and the output is y(t). Figure 1 - 1 5b 
shows the equivalent gain representation of Q^. The equivalent gain 
of Qp is defined as 



(z) = 


V z) 

t<x (z) 


(1-51) 


where <I> ( z ) denotes the z-transform of the autocorrelation function of 

x(t), and $ (z) is the z-transform of the crosscorrelation function of 

xy 
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x(t) and y(t). In order to derive K^z), we consider that the input 
to Qp is a stationary Gaussian process with zero mean so that the 
probability function is 


P(x) 


exp(-x 2 /2 0 2 ) 
/tn a 

x 


where a is the standard deviation of x(t). 

X 

The crosscorrelation function of x(t) and y(t) is given by 



(x) 


°° 

Q D (A)p(X)pXdX 

-00 


where 


Q D Cx(t)] = y(t) 


is the mathematical description of the quantizer, and 


*xx ( t) 
P 2 


Thus, 


f 00 

$ X y( T ) = ~2 ^ xx (t) j %0)p(A)AdA 


Taking the z-transform on both sides of Eq. (1-56) and rearranging, we 


* xv (z) 1 


Q D (X)p(X)XdX 


(1-52) 


(1-53) 


(1-54) 


(1-55) 


0-56) 

have 

(1-57) 


Although the equivalent gain ^(z) is indicated as a function of z, 
the nature of the right-hand side of Eq. (1-57) implies that it is 
always a constant. It is apparent that although our interest is centered 
on the quantizer nonlinearity, in general, the definition of ^(z) can 
be applied to any comman nonlinearities found in control systems. 

With reference to the quantizer characteristics of Fig. 1-7, 


y(t) = 

0 

-h/2 < x < 

h/2 


y(t) = 

-h 

-3h/2 < x < 

: -h/2 


ri- 

ll 

h 

h/2 < x < 

3h/2 


y(t) = 

- (N-l )h 

2 ( N- 1 ) h . 
2 

X £ - 

(2N-3)h 

2 

y(t) = 

(N-l)h 

2 ( N- 3 ) h , 
2 1 

x<Mi 

y(t) = 

-Nh 

-C» < 

X £ - 

(2N-1 )h 

2 

y(t) = 

Nh 

2 ( N- 1 ) h . 
2 - 

X < oo 



where N is a positive integer. 
Equation (1-57) gives 




r- ( 2N- 1 )h/2 

i (-Nh)p(A)AdA 


+ 


+ 


r (2N-3)h/2 

- ( N- 1 )hp(A)AdA + ••• + 

- (2N-1 )h/2 

f 3h/2 f (2N-l)h/2 

hp( A ) AdA + • • • + 

J h/2 ( 2N-3 )h/2 


-h/2 

-hp(A)AdA 

J-3h/2 


(N-l ) hp ( A ) AdA 


+ 


Nhp(A)AdA 

(2N-1 )h/2 


( 1 - 
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Or, 


K eq (z) * -7 

a 

x 


+ 2Nh 


2h 


3 h/2 
h/2 


p(X)XdX + 


2(N-1 )h j 


(2N-1 ) h/2 
( 2N-3 )h/2 


p(X)XdX 


(2N-1 ) h/2 


p(X)XdX 


0-59) 


Since all the integrals in the last equation have the same integrand, 
we evaluate one of the integrals as follows: 


r 3h/2 

h/2 


p(X)XdX = 


f3h/2 exp(-X 2 /2a 2 ) 


h/2 


-= — XdX 

/2TT a„ 


(1-60) 


Let 


Then 


u = X 2 /2a 2 


du = 


XdX 


or 


XdX = a du 

A 


and 


3h/2 

h/2 


p(X)AdX = 


1 


2 2 
9Ji /8a 


? * exp(-u)cr du 

/2ir J h /8a 2 X 


/2? a x 


2 2 } 
exp(- ) + exp(- ^) 


8a 


8a? 


J 


d-61) 


Thus, the equivalent gain for the quantizer with 2N levels of quantization 
i s 
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K eq (z) "/t r a x 


exp(- + exp(- ^4-) + exp(- ^4-) + ••• + exp(- 
8 a 8a 8a z 

XXX 


(2N-T) 2 h 2 
8a 2 


x 

(1-62)' 


Although the equivalent gain has been derived from Q^, the 
response of the closed-loop system of Fig. 1-14 cannot be determined 
by using K (z) directly, since a , the standard deviation, of x(t) 

GCj X 

is not known. The analysis procedure is outlined as follows: 

For a given a c , which is the rms value of the input, 

(1) The equivalent gain K (z) is computed using Eq. (1-62) for 

C 

various values of a . 

A 

(2) The values of a that correspond to the various values of K (z) 

x eq 

obtained in step (1) are calculated from 


a x = 


1 


2?r j 


-1 


-1 


•- 0 M (z)M (z" )4> (z)z" dz 


cc 


1/2 


(1-63) 


where 


M x (2) 


1 

1 + K eq (z)G A (z)G B (z) + G c (z) 


(1-64) 


and 


’■cc* 2 * 


= a. 


It is assumed that the input of the system is a white noise. 

The solution of o for the given a is determined when there is 
X c 

a match between a from Eq. (1-63) and that used in step (1). 

A 

(3) The rms attitude response of the system is determined from 
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On = 


2irj 


(|) M(z)M(z _1 )$ cc (z)z _1 dz 


1/2 


(1-65) 


using the k^z) which corresponds to the a x obtained in 
step (2). 

For the LST system, a quantization level of h =0.003 is selected 

for Q d> The quantizer is assumed to saturate after 5000 increments; that 

is, N = 5000 in Eq. (1-62). A digital computer program is prepared 

which automatically cuts off the series of Eq. (1-62) when an additional 

term is contributing less than 10'^ to the entire result. The sampling 

period T is chosen to be 25 msec, and the system transfer functions G^(z), 

Gg(z), and G c (z) are given by Eqs. (1-38), (1-36), and (1-37), respectively. 

Following the procedure outlined above, the results of the analysis with 

several values of a are tabulated in Table 1-4. 

c 

The results in Table 1-4 show that for a given set of quantization 
and saturation levels, both small and large input signals cause the 
quantizer to act as an attenuator. The true characteristics of the system 
as a function of the input are displ ayed by normalizing Og. The last 
column in Table 1-4 represents the normalized output Og^, which is 
defined as 


°BN 


g B 

a Bl a c 


where 


( 1 - 66 ) 


a Bl a B at a c ‘ 1 -• 


0-67) 
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Table 1-4. 



Input 


No. of terms 

Output 

Normal i zed 

0 

a 

u used in K (z) 

K eq 

O Q 

output 

c 

X 

eq ^ 

B 

°BN 


5x 10‘ 4 

5.318 x 10“ 4 

0.075 

2 

5.20 x 10~ 5 

0.4 

10“ 3 

1.086 x 10" 3 

0.797 

3 

2.40 x 10 -4 

0.92 

10 -2 

1 .093 x 10" 2 

1.0 

19 

2.68 x 10" 3 

1 .0 

10" 1 

1.093 xl0 -1 

1.0 

176 

2.60 x 10~ 2 

1 .0 

1 

1.093 

1.0 

1753 

2.60x 10 _1 

1 .0 

10 

10.87 

0.836 

5000 

2.383 

0.916 

50 

53.37 

0.22 

5000 

6.88 

0.529 


Figure 1-16 shows the plot of ^(z) versus a^. The plot has the 
significance of a "statistical describing function" of the quantizer. This 
plot shows that for small inputs, when the input magnitude is comparable to 
the quantization level, h, the gain drops below unity, and the quantizer is 
attenuating the signal. For larger inputs the quantizer appears as a unity- 
gain element in a statistical sense, and for very large inputs, where the 
input magnitude is comparable to maximum output (N levels of h), the 
quantizer gain again reduces, and the input is attenuated. 

Figure 1-17 gives a plot of a versus a which represents the solution 

X c 

of Eq. (1-63) and step (1) for the range of ct c considered. 

Figure 1-18 shows a plot of 0 ^ versus o . It is seen that the normalized 
output reduces from unity in the ranges of when ^(z) is less than 1. 

As an interesting comparison, the statistical method of Section 1-5 
assumes ^(z) = 1 for all values of o Q and, consequently, it always yields 























46 0780 



33. 3h 333h____ 3333h_J~ 33333h- 





















41 


a result of Og^ = 1. In'this sense that method is more conservative and does 
not recognize the severity of the nonlinearity due to the quantizer at 
very small and very large input magnitudes. 
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2. Stability of the Analog Low-Cost Large Space Telescope Due to 
Quantization 

2-1 . Introduction 

The objective of this chapter is to conduct an investigation on 
the pointing stability of the low-cost Large Space Telescope (LST) 

System under the influence of quantization at various locations of 
the system. Only the analog model of the LST system is considered 
in this chapter. The stability of the digital LST system with quan- 
tization is considered in Chapters 5 and 6. 

The block diagram of the simplified analog LST system with 
quantizers is shown in Fig. 2-1. 

The describing function method is used to determine the condition 
of self-sustained oscillation in the LST system due to the effect of 
each of the three quantizers acting alone. The present analysis 
considers only one quantizer at a time. 

Let the analog describing functions of the quantizers Q^, Q^., and 
Q R be represented by N Q , N-p and N R , respectively. In general, the 
describing function of a quantizer nonlinearity is a function of the input 
amplitude, E; quantization level, h; and the number of quantization 
levels, n, which depends on E. 

For the LST system shown in Fig. 2-1, the "characteristic equations" 
of the system when each one of the quantizers is acting, are given as 
follows: 
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Displacement Quantizer Q^: 


V 3 + n d (k p s + K I ) + k r $ 2 = 0 
Torque Quantizer Q^: 


J v s 3 + N t (K r s 2 + Kp's + Kj) = 0 


(2-1) 


( 2 - 2 ) 


Rate Quantizer Q R : 

J v s 3 + K p s + Kj + N r K r s 2 = 0 (2-3) 

These equations can be conditioned by dividing both sides of the 
equations by the terms that do not contain the describing function, so 
that the stability equation is expressed in the form of 


1 + NG(s) =0 (2-4) 

where N is the describing function, and G(s) is a linear transfer 
function. The condition of self-sustained oscillations is found by 
investigating the possible intersections between the trajectories of 
-1/N and G(s) in the complex plane. 

For the equations given in Eqs. (2-1), (2-2) and (2-3), the equiva- 
lent transfer functions are 


Gp(s) 


V 1 K i 

s 2 (0 v s ♦ K r ) 


(2-5) 
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K r s + K s + K t 

g t (s) = P L 

1 J s J 

V 


G r (s) = — — -S 

V + K P S + K I 
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2-2. Analog Describing Function of the Quantizer Nonlinearity 

Consider that a quantizer has the input-output relation as shown 
in Fig. 2-2. Let the input to the quantizer be a sine wave, 


e(t) = E sin ait 


(2-8) 


The quantization level is h, and let the magnitude of E be such that 


(2n - 1 )h £ (2n + 1 )h 

2 - 2 


(2-9) 


where n is a positive integer. A typical output of the quantizer is 
shown in Fig. 2-3. 

The fundamental component of the Fourier series expansion of y(t) 
is 


y-j (t) = Y 1 sin ait 


( 2 - 10 ) 


where 


4 

7r 


tt/2 


y(t) sin ait doit 


( 2 - 11 ) 


Evaluating the last integral, we have 




a 4 ffr/2 

y(t) sin ait doit + - 
a l 


nh si n ait doit 


a 


n 


4 

” i=i 


a. 


l+l . . . . , 4 

ih sin ait doit + - 

71 

“i 


77/2 


nh sin ait doit 


ot„ 


( 2 - 12 ) 
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Then, 


n 1 

- — y i[cosa.,-, - cos a-] + ^-cosa 
7T . L , 1+1 1 TT n 

1=1 


Si nee 


sina^ = 


h 

2E 


(2i - l)h 
sina^ = A ^ 


cos a-j = /l - 


2E 


cos a . = / 1 - 


f(2LrJLLh 


i2 


2E 


Expanding the right-hand side of Eq. (2-13), we can show that 


4h 


n 


Y, = — 7 cos a. 
1 * i=i 


Substituting Eq. (2-17) into Eq. (2-18), we have 


Yn = — l (l - 
1 " i=i 1 


f(2j - 1 )hl 


2E 


J J 


1/2 


The describing function of the quantizer is written as 


f 1 - 


i=l 


2E 


2 ) 1/2 


(2-13) 

(2-14) 

(2-15) 

(2-16) 

(2-17) 

(2-18) 
■ (2-19) 


(2-20) 
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Figure 2-4 shows the plot of N(E/h) as a function of E/h. We 
shall show in the following that 


tun N ( E/h ) = 1 
E 


Letting x = E/h, Eg. (2-20) is written 


M(X) .l =] V ~ 


f2i - ll 

2 X 


2 ) 1/2 


(2-21) 


( 2 - 22 ) 


Expanding the quantity inside the summation sign in the last 
equation, we have 


N( X ) 


~ l 

i i 1 


Or 


1 1 

(i j_i 

' i| 

n 

2 

.X " 2 Xj 

' 8| 

X 1 

1 

roj 

xj 


3 

fl J .1 

6 15 | 

[i JJ 

8 

48 

X 1 
1 

rol 
X | 

384 J 

.X ' 2x, 

> 


(2-23) 


N( x>=^ 




15 


384 x c 




■ (2-24) 


Taking limit as x approaches infinity and n approaches infinity, in 
the last equation, we get 
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tun N ( x ) s tin — l 
nx i=l 


X^ 

n-x« 


X^° 

n-*» 


.2 


.4 


1 - 


3i 


.6 


1 5i 


8 


2 X 8x 48x 384 x 


8 


tun — 

•nx 

X -K» 

n-»°° 


n - 


Y_ H_ 
2 3 


2x 


1 n 
4 5 


15 n- 


8x 


48 X 


6 7 


384 x 


8 9 


4 

TT 


1 - 1-1 
1 2 3 


11 A. 1 

8 ’ 5 " 48 ‘7 


15 1 


384 9 


4 

TT 


1 A 2 , 

/ 1 - x ax 


4 

TT 


/ 2 
X xl - X 


+ sin -1 x 


=1 (2-25) 

Since the describing function of the quantizer nonlinearity is 
always a real number, the function -1/N(E/h) in the magnitude (db) 
versus phase coordinates will lie on the -180-degree axis for all 
values of E/h. The plot in Fig. 2-4 shows that the magnitude of 
-l/N(E/h) is infinite for 0 £ E/h < 0.5. For 0.5 £ E/h £ 0.707, the 
plot of -l/N(E/h) extends from infinity to -2.09 db along the -180-degree 
axis. Over the range of 0.707 £ E/h £ 1.5, -l/N(E/h) extends from 
-2.09 db to 1.9 db along the -180-degree axis, etc. As E/h approaches 
infinity, the plot of -l/N(E/h) is reduced to the zero-db point. 

Figure 2-5 shows the plot of N(E/h) in magnitude versus E/h and shows 
the multivalued property of the function. 
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N 

0 
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1.5 


1.25 -1 


1.7 5 

1.0 H 


U 0.53 


h 0.55 


- E/h = 0.6 

- 0.65 

E/h = 0.707 


-180° 

PHASE 

. Analog describing function for the quantizer nonlinearity. 
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2-3. Self- Sustained Oscillations of the Analog LST System with 
Quantization 

The transfer functions of the three quantizers, Q^, Qy, and Q R , 
given by Eqs. (2-5), (2-6), and (2-7), are plotted in Fig. 2-6, together 
with the -l/N(E/h) trajectory. The following system parameters are used: 

Kp = 1 .65 x 10 6 K r = 3.71 x 10 5 

Kj = 7.33 x 10 5 J y = 41822 

Figure 2-6 shows that since the plot of Gp(s) does not intersect 
the -1/N(E/h) trajectory in the finite domain, the quantizer Q D will 
not cause any self-sustained oscillations in the LST system. 

Both the curves for Gy(s) and G R (s) intersect the -180-degree axis 
and the -l/N(E/h) trajectory at 25 db. 

Thus, 

20 £og 10 |l/N(E/h)| = 25 db (2-26) 

which gives 

N ( E/h ) = 0.05012 •' (2-27) 

Figure 2-4 shows that for this value of N(E/h), 

£=0.5 (2-28) 

or 


E s o.5 h 


(2-29) 


DECIBELS 








This means that the quantizers and Q R each independently may 
cause self-sustained oscillations to occur in the LST system. Further- 
more, the self-sustained oscillations will have an amplitude equal to 
one-half the level of quantization h. Another point of interest is 
that the intersect between G-j.(s), G R (s), and -l/N(E/h) corresponds to 
n = 1 , so that the quantizers Q^. and Q R are essentially acting as a 
simple relay with dead zone. 

Figure 2-6 shows that the frequency of the oscillations caused by 
Q-j. is 1.5 rad/sec, and that caused by Q R is 6.3 rad/sec. 

It is interesting to point out that the analysis in Chapter 1 shows 
that the displacement quantizer Q q produces far greater quantization 
error than the torque and rate quantizers. However, the describing 
function analysis shows that Q D does not cause self-sustained oscillations, 
whereas and Q R may excite oscillations with amplitudes equal to one- 
half the level of quantization, h/2. All these factors must be taken 
into consideration when selecting the quantization level. 
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3. Stability of the Analog Low-Cost Large Space Telescope with 
Reaction Wheel Friction Nonlinearity 

3-1 . Introduction 

The objective of this chapter is to conduct an investigation on 
the pointing stability of the analog model of the low-cost Large Space 
Telescope (LST) System with the reaction wheel frictional nonlinearity. 

Although the LST system is digital, the stability study on the 
analog model will establish a limiting case when the sampling period 
becomes very small, thus providing a check on the results of the digital 
system. 

The block diagram of the digital low-cost LST system, including the 
reaction wheel dynamics, is shown in Fig. 1-1. It has been shown in 
Section 1-1 that the dynamics of the LST system can be simplified. 

Figure 3-1 shows the block diagram of the simplified low-cost LST 
system with the reaction wheel nonlinearity. 

The A of the system shown in Fig. 3-1 is 

A = 1 + K r G 4 + G^Gg + NG 6 G 7 (3-1) 


where N denotes the continuous-data describing function of the reaction 
wheel frictional nonlinearity. 

Substitution of the expressions of , G^, G^, Gg, and Gy into 
Eq. (3-1 ) , we have 


K R K I 1 

A=l+^r+(K p +-f)-^ 


J s 
v 


J s 

V 


1 

J s 2 
J RW S 


(3-2) 
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3-2. Condition of Self-Sustained Oscillations in the Analog 1ST System 
with Reaction Wheel 


It has been shown [1] that the nonlinear frictional characteristics 
of a reaction wheel can be described by the Dahl model. Therefore, the 
analog describing function of the CMG frictional nonlinearity derived 
in [5] can be directly utilized. 

The equivalent transfer function that N sees is determined from 
Eq. (3-3) by dividing both sides of the equation by the terms that do 
not contain N. We have 


NJ s 

1+ = ^ =0 (3-4) 

J RW (J v S + K R s + V + K I> ’ 

Thus, 

J s 

S M (s) = 3 ~2 < 3 - 5 > 

q J RW<V + K R s + V + K I> 

The condition of self-sustained oscillation is investigated by 
plotting G e q(jto) and -1/N in the magnitude (db) versus phase coordinates. 
Figure 3-2 shows the plot of G^jcj) of the LST system with 

K r = 3.71 x io 5 

K = 1.65 x 10 6 
P 

Kj = 7.33 x 10 5 
J v = 41822 



and the plots of -1/N with 


(i) Y = 8470 (ii) Y = 84700 ( i i i ) y = 

t gfo - °- 0424 t gfo = °' 0424 t gfo = 

Since the G ec| ( joo) plot never enters the region bounded by -180° 
and -270° in which the plot of -1/N lies, the analog LST system will 
not have self-sustained oscillations, due to the reaction wheel 
frictional nonlinearity. 


8470 

0.424 
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4. Stability of the Digital Low-Cost Large Space Telescope With 
Reaction Wheel Friction Nonlinearity 


4-1 . Introduction 

In this chapter we shall investigate the condition of self- 
sustained oscillations of the digital low-cost LST system with reaction 
wheel friction nonlinearity. 

The block diagram of the simplified LST system is shown in Fig. 

4-1. The equations written for the outputs of the samplers are: 

$ e (z) = -G A (zH e (z) + N(z)G b (z)0 rw (z) (4-1) 

0 RW (z) = G c (z)d> e (z) - N(z)G d (z)0 rw (z) (4-2) 


where N(z) denotes the discrete describing function of the reaction 
wheel friction nonlinearity. Since the friction nonlinearity of the 
reaction wheel can be represented by the Dahl model, N(z) is identical 
to the discrete describing function of the CMG nonlinearity [5]. 

The transfer functions in Eqs. (4-1) and (4-2) are 


G A (Z) = ft 


G ho G l G 2 G 3 
1 + K R G 2 


G B (z) = ft 


G ho G 2 G 3 
1 + K r G 2 


G c (z) vft 


EhoW 

1 + K R G 2 


(4-3) 


(4-4) 


(4-5) 
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g d ( z > = 2r 


fG ho G 4 G 5 


1 + K r G 2 


The signal flow graph representing Eqs. (4-1) and (4-2) is shown 
in Fig. 4-2. The A of the system is obtained from Fig. 4-2, 


A - 1 + G^( z) + N(z)[G^(z)Gp(z) - Gg(z) G c (z) + Gp(z)] 


Setting A to zero, the last equation can be written in the form 
of 


1 + N(z)G eq (z) = 0 


where 


G (z) 
eq ' 


G/\(z)Gq(z) - G b (z)G c (z) + G d (z) 


1 + G a (z) 


The individual transfer functions are evaluated as follows: 


g a ( z > ■ <’ - z_1 )f y 

v w 




s 2 (s + a) s 3 (s + a) 


K p G B (z) + J 


I 


f T 2 (z + 1) 


J + 1 


z - 1 


[2a(z - l) 2 a 2 (z - 1) a 3 a 3 (z - e -aT )j 


a - K d /J 
R v 


~ K n 


1 - e 


-aT 


Z ' 1 a(z - e’ aT ) 


(4-6) 


(4-7) 


(4-8) 


(4-9) 


(4-10) 


(4-11) 
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u) . /- V zl 


G u) = /- 


(4-UV 

( 4 - 13 ) 
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4-2. Self- Sustained Oscillations in the Digital LST System With 
Reaction Wheel Nonlinearity 

Figure 4-3 shows the plots of. G (z) i n Eq- (4-9) for various 
values of n with T as a parameter. The integer n and sampling period T 
are related to the frequency of oscillation by the following equation: 


to 

c 


2tt 

nT 


The following system parameters are used: 


J 


K 


P 


RW. 


= 1.65 x 10 6 
= 7.33 x 10 5 
= 0.2 


K r = 3.71 x 10 5 

J =41822 
v 


(4-14) 


In Fig. 4-3 the curve for n = 2 extends up to approximately 30 db 
at T = 0.7 sec. Also, as n is increased the curves for G ec) (z) approach 
the curve for G e ^(s) presented in Fig. 3-2. The plots of G eq ( z ) 
together with the plots of -1/N(z) allow the study of self-sustained 
oscillations in the digital LST system. 

Figures 4-4 through 4-6 show the plots for -1/N(z) for various n. 
In these plots y = 84700 and T^q = 0.424. In all cases, the magnitude 
of the lowest point of the -1/N(z) curve as E approaches 0 is given by 


tan 

E-MD 


1 

wr 



(4-15) 


With the given parameters this point is at -83.66 db. The curves 
for n = 2 are shown in Fig. 4-4 and consist of two straight lines at 
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-180° (for 0 £ 4) < tt/2 ) and -360° (for tt/2 < <J> £ it). The plots for 
n = 3 and n = 4 are shown in Figs. 4-5 and 4-6, respectively. Several 
values of <p are plotted in each case to illustrate the effect of the 
phase of the input signal. It should be noted that for odd n the curves 

repeat every 180/n degrees starting from <f> = 0, and for even n the curves 

repeat every 360/n degrees starting from <j> = 0. As the input amplitude E 

goes to infinity, the curves also go to infinity and span a region 180/n 
degrees or 360/n degrees wide for odd or even n, respectively. This 
region is centered about the -270° line. Thus, as n goes to infinity 
these curves approach the -1/N curve of the continuous system shown in 
Fig. .3-2. 


For stability analysis it is sufficient to consider only the bounds 
of the - 1 / N ( z ) plot for a fixed n. Self-sustained oscillations can 
occur if the G(z) curve corresponding to the same n intersects with 
the F(z) plot. Figures 4-3 through 4-6 show that self-sustained oscillations 
can readily occur in this system with the choice of T determining the 
possible n values which can exist. As the critical regions for higher 
values than n = 4 can be easily visualized, it is apparent that this 
system will have self-sustained oscillations for integral values beyond 
n = 4. However, as n increases, theamplitude of oscillation decreases 
and the oscillation will eventually cease as T gets smaller and smaller 
and n gets larger and larger. 

If y and Tg^Q are increased then the curves of -1/N(z) move up with 
the end point shifting according to Eq. (4-15). Figures 4-7 through 4-9 
show the plots of -1/N(z) for n = 2, 3 and 4, respectively, with y = 8470 
and Tgpg = 0.0424. In this case the lowest point is at -23.66 db and 




Figure 4-3. G(z) curves of the low-cost LST with frictional nonlinearity. 













































Figure 4-7. Discrete describing function for the frictional nonlinearity 
n = 2. J * 



































self-sustained oscillations are possible only for n = 2 over a fixed 
small range of sampling periods. If T is chosen to exclude the critical 
region, self-sustained oscillations can be avoided. 
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5. Discrete Describing Function of A Quantizer 
5-1 . Introduction 

The effects of quantization and the study of quantization error 
conducted in Chapter 1 are all based on the assumption that the digital 
system with quantization is stable and free from sustained oscillations. 
Since a quantizer is a nonlinear element, it can cause self-sustained 
oscillations. In Chapter 2 the condition of self-sustained oscillations 
in the analog low-cost LST system with quantizer is studied by use of 
the continuous-data describing function. However, in reality, the low- 
cost LST system is a digital system. The interaction between the sampling 
operation and the quantizer will bring about phenomena which can be 
grossly different from that in an analog system. Therefore, it is 
essential that the discrete describing function (DDF) of a quantizer be 
derived. To the authors' knowledge the DDF of a quantizer has not been 
derived in the past. 

Figure 5-1 shows the input-output characteristics of a quantizer. 

The input of the quantizer is denoted by e*(t), and the output by y*(t). 

It is assumed that the input of the quantizer is the output of an ideal 
sampler. Therefore, e*(t) and y*(t) are trains of impulses. Furthermore, 
it is assumed that the input to the sampler, e(t), is a cosine wave, and 
thus the amplitude of e*(t) is modulated by a cosine wave; that is, 

e{t) = E cos (wt + <|>) (5-1 ) 

co 

e*(t) = E l cos (koaT + <f>)6(t - kT) 
k=o 


(5-2) 
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where E is a constant, oj is the frequency in radians per second, T is 
the sampling period in seconds, and 6(t) denotes the unit impulse 
function. The z-transform of e*(t) is 


E(z) 


Ez 


z - 2z cos aiT + 1 


[(z - cos coT) cos 4> - sin wT sin <j>] 


(5-3) 


the z-transform of y*(t) is denoted by Y(z). 

The discrete describing function of the quantizer nonlinearity is 
defined as 


N(z)=|||| (5-4) 

As in the case of the relay-type nonlinearity, we assume that 
because of the periodic nature of the sampler, e*(t) and y*(t) are all 
periodic functions of period nT, where n is a positive integer greater 
than or equal to two. Thus, 



(5-5) 


or 

^ = ~ (5-6) 

n = 2, 3, 4 


The DDF N(z) is incorporated in the "characteristic equation" of 
the system. 


1 + N(z)G eq (z) = 0 


(5-7) 


for the determination of the condition of self-sustained oscillations, 
where G (z) denotes the linear transfer function which the quantizer 
nonlinearity sees. Graphically, the condition of self-sustained 
oscillation characterized by the period nT is determined by the inter- 
sections of the ^(z) trajectories with the critical regions of -1/N(z), 
all for the same n. Therefore, the DDF problem involves the determination 
of the critical regions of -1/N(z) for the quantizer for n = 2, 3, ... 


81 


5-2. The DDF of A Quantizer for n = 2 


In order to illustrate the derivation of the DDF of the quantizer 
nonlinearity, we shall first consider the case of n = 2; that is, the 
self-sustained oscillation is characterized by the period which is equal 
to twice the sampling period. 

For this mode of oscillation, the waveforms of e(t), e*(t), and 
y*(t) are shown in Fig. 5-2. 

For n = 2, Eq. (5-6) gives ojT = tt; Eq. (5-3) becomes 

E(z) = E - Z - C Q^ ' (5-8) 

Z T I 

The z-transform of y*(t) is written 

YU) =VtT (5-9) 


where k is a positive integer. In this case it is assumed that the 
value of E is constrained by the following equation: 

i2L^ <Ecos+< il!LiJih 

Using Eqs. (5-8) and (5-9), we have 


1 E cos <(> 

NUT ‘ kh 


(5-11) 


For a given set of values of k, h, and <}>, the constraints on the 
values of E for n = 2 are 
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r = (2k + jjh 

max 2 cos 4> 

c = (2k - 1)h 

min 2 cos 4> 

Substituting E max and E^ into Eq. (5-11), we have 


Nib 


(k) 


max 


NliT 


(k) = - 


min 


(2k + 1) 
2k 


(2k - 1) 
2k 


(5-12) 

(5-13) 


(5-14) 


(5-15) 


The last two equations indicate that the critical region for -1/N(z) 
for n = 2 is the line which extends from -(2k - 1 )/2k to -(2k + 1 )/2k 
on the negative real axis in the polar coordinates. 

It can be shown that 


N(z) 


( k+1 ) 

max 


NTz) 


(k) 

max 


Since 


1 

Niry 


(k+l ) 
max 


2(k + 1 ) + 1 
2(k + 1) 


+ 


1 

2(k + 1 ) 


and 


N(z) 


(k) 

max 


2k + 1 
2k 


1 + W 


(5-16) 


(5-17) 


(5-18) 


thus, Eq. (5-16) is verified. 
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Similarly, using Eq. (5-15) we can show that, in general, 


NuT 


(k+1) 


min 


NUT 


(k) 


nn n 


Thus, let S(k) be the set which is bounded by 


Wz) 


(k) 


and 


max 


N(zT 


(k) 


min 


and S(k+1) be the set which is bounded by 


- nTzT 


(k+1) and 


max 


N(TT 


(k+1) 


mm 


then 

S ( k+1 ) C S(k) 

for k=l,2,3, ... . 

Equations (5-14) and (5-15) also show that 


(5-19) 


(5-20) 


tim 

k-»o 


T 


(k) 

max 


= -tim 
k-K) 


NlzT 


(k) 


min 


= -1 


(5-21) 


This result implies that as the number of quantization levels increases, 
the quantizer characteristics approach a linear gain, and the critical 
region reduces to the (-1, jO) point in the complex plane. 

The critical region for n = 2 is shown in Fig. 5-3. 
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5-3. The DDF of A Quantizer for n = 4 

Before embarking on the derivations of general expressions of 
-1/N(z) for all n > 2, the case for n = 4 will be considered as an 
illustrative example. 

For n = 4, the input pulse train of the quantizer can assume a 
maximum of two different pulse amplitudes, k Q and k^ . The mode of 
oscillation is thus characterized by A = (k Q , k-j ) , where k Q and k-j are 
positive integers. Figure 5-4 illustrates the input and the output of 
the sampler when the former is shifted through 360 degrees. It is 
observed from these waveforms, as well as will be verified later by 
equations, that the critical region for n = 4 repeats every 90 degrees 
for <}), so that only the range for -45° </<(>< 45° needs be considered. 

For n = 4, cos wT = cos 90° = 0, z = j. Thus, Eq. (5-3) becomes 


E(z) = — — (z cos <p - s i n ) (5-22) 

z c + 1 

With reference to Fig. 5-4, the expressions of the input to the 


quantizer for the various 


-90° £ 4> £ 0° 


-180° < 4> £ -90° 


0° < cp < 90° 


ranges of cp are written 
h(k z + k,)z 

V(2) = —2j !— 

z + 1 

h(-k z + k,)z 

Y(z) = % — 

z~ + 1 

h(k z - k,)z 

V(z) = ^ 5— 

z + 1 


follows: 

(5-23) 

(5-24) 


(5-25) 




2T 


3T 


4T 



Figure 5-4. Input and output signal wavef 
the quantizer as <f> varies thn 



t 


0° < 4> < 45° (k ] < k 0 ) 



t 


* = 45° (k ] = k Q ) 



t 


45° < <t> < 90° (k > k ) 
1 o' 



t 


<P = 90° (k Q = 0) 



t 90° < ij> < 135° (k Q < k 1 ) 



♦ = 135° (k o = k^ 



t 135° < $ < 180° (k Q > k ] ) 


7TT 


t 


<p = 180° (k 1 = 0) 


of sampler preceding 
360°. 
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h{-k n z - k,)z 

90° < <j> < 180° Y(z) = % 1 — 

+ 1 


(5-26) 


The constraints on the magnitude of E over these ranges of <j> are 
tabulated as follows: 



E cos <p | 


(2k o + l)h 
_ 


and 


( 2 k-j - l)h 
2 


(2k, + 1 )h 
|Esin<j>| £ p 


(5-27) 


(5-28) 


These inequality conditions lead to maximum and minimum bounds on E. 

When these E „ and E . are substituted into 
max min 


1 

N(z) 


E(z 

Y(z 


(5-29) 


and along with Eqs. (5-23) through (5-26), the equations that define 
the boundaries of the critical regions are obtained as 


NlzT 


and 


max 


NUT 


min 


For instance, for -90° £ <j> £ 0°, from Eq. (5-27), 


E . 


(2k - l)h 


min 1 2 cos <(> 


(5-30) 


Using Eqs. (5-22), (5-23), and (5-29), we have 
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NUT 


2k o- 1 


= 2tk r+ k ) (tan * _ 

min 1 K 1 1 


(5-31) 


Rationalizing the last equation, we have 


NTYT 


2k o " 1 

: o — [ ( k, tan 4) - k ) + j(-k tan 4> 

mini 2(k£ + k 2 ) 1 0 0 


k-j )] 


(5-32) 


Let 


m 


min 1 


= R + jl 


(5-33) 


where 


2k - ! 

R - o p (k-i tan 4> - k ) 
2(k 2 + k 2 ) 1 0 


(5-34) 


2k - 1 

I = - — y 2~ (-k tan 4 ) - k, ) 

2(k ‘ + kf) 0 1 


(5-35) 


Then, solving for tan 4) from Eq. (5-34) gives 


R + k C 

tan< » = -Tqc“ 


(5-36) 


where 


C = 


2k o " 1 
+ k (> 


(5-37) 


90 


Substitution of Eq. (5-36) into Eq. (5-35) yields after simplification, 


k 2k - 1 

r - _ 0 D _ 2 

1 k K 2k 

I ^ K 1 


(5-38) 


which is the equation of a straight line in the complex -1/N(z) domain. 
Similarly, from Eq. (5-27), 


( 2 k o + ! )h 

"max 1 2 cos <j> 


-90° < (j> < 0 C 


(5-39) 


Following through the same procedure as described above, we have 


N(z) 


2k + 1 


max 


, ° 2(k o j V T7! < tan f - 


(5-40) 


which is represented by a straight line with the equation 


i- !!kll 

k i 2k i 


(5-41) 


The boundaries for -1/N(z) when the constraint equation of Eq. (5-28) 
is used are denoted as 


1 


N(z) 


and 


1 


max 2 


nTzT 


min 2 


The equations which define the boundaries of - 1 / N ( z ) for the entire 
range of <j> are tabulated below: 


■90° < <p < 0° 
90° < <|> < 180 c 


k 2k - 1 

T - _ _2 D _ 0 

1 k ] K 2k 1 


(min 1 ) 
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k 2k + 1 

T = _ _2. D _ 0 

k l 2k l 


(max 1 ) 


k, 2k, - 1 
0 0 


(min 2) 


k, 2k + 1 

I= ir R+ 4k— 

0 0 


(max 2) (5-42) 


-180° £ 4> £ -90 c 
0° < 4) < 90° 


k 2k - 1 

T - 0 D J. 0 

’ 2k l 


(min 1 ) 


k 2k + 1 

I = ^ R + -V 


(max 1 ) 


k, 2k, - 1 
o o 


(min 2) 


k, 2k + 1 

1 ■ - IT * - 4k— 

o o 


(max 2) (5-43) 


These equations show that for the oscillatory mode of n = 4 and 
A = (k Q , k-j), the critical regions are bounded by straight lines, and 
that it is sufficient to consider the range of -90° £ 4> £ 90°, 

However, Fig. 5-4 shows that k Q £ k-j for -45° £ 4> £ 45° and k Q £ k-j 
for 45° £ |4>| £90°, thus, the critical regions for modes A = (k Q , k^ ) 
and A = (k-j, k Q ) are identical. It is necessary only to consider the 
range of <p from -45° to +45°. In general, one does not have to use the 
constraints on 4>, as the intersects of the eight equations in Eqs. (5-42) 
and (5-43) will naturally define the critical region. 
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Modes A = [k o> 0) and A = ( 0 , ) 

Figure 5-4 shows that when <}> = 0° or 180°, k-| = 0, and the mode is 
described as | = (k Q , 0). Similarly, when (j> = ±90°, the mode is 
A = (0, k 1 ). 

For the mode A = (k Q , 0), Eqs. (5-42) and (5-43) are reduced to 
2k - 1 

R = gk ( min 1 ) 

P 


2k i + 1 

R = ~r (max 1 ) 


(min 2) 


(max 2) (5-44) 

These four lines define a square as shown in Fig. 5-5. 

i 

For the mode A = (0, k-j), it is simple to show that the critical 
region is described by equations of the same form of Eq. (5-44) with 
k Q replaced by k-j . Therefore, the square region shown in Fig. 5-5 
is also for the liiode A = (0, k^ ) with k Q replaced by k-j . 

Mode A - [k Q , fej ) , k Q = fej 

When <j> = ±45°, ±135°, k Q = k^, as shown in Fig. 5-4. The critical 
region is now described by the following four equations: 

2k - 1 
o 

2k o 



I = -R - 
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I 


-R - 


2k 0 + 1 
2k o 


I = R 


+ 


2k o 


I = R 


+ 


2k + 1 
o 

2k 

o 


l 

The critical region is again a square, as shown in Fig. 5-5. 

In general 1 , the critical region for the mode A = (k Q , k^ ) is . 
defined by the intersects of the eight equations in Eqs. (5-42) and 
(5-43). As an illustration, Fig. 5-6 shows the critical region for 
A = (2, 1) which is also for A = (1, 2). Notice that this critical 
region is bounded by that of A = (2, 0) or A = (0, 2). 




i 


4 



(5-45) 
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5-4. A General DDF of a Quantizer for n > 2 

In this chapter we will generalize the results of the previous 
sections and derive a DDF of a quantizer which is valid for all n > 2. 
Thus, for an oscillation of any order, n, we seek expressions which 
will define the critical regions of -1/N(z) for all possible modes 

i 

that can exist. 

The derivations for the case of odd n and for the case of even 
n differ substantially and are therefore considered separately. 

i 

A. Even n 

I 

When n is even, the input pulse train to the quantizer can 
assume a maximum of n/2 different amplitudes. The mode of oscillation 
is thus characterized by 


A " (k o’ k T k 2 ’ k n 1 -l ) 

i'; 


where 

n-j = n/2 


(5-46) 


(5-47) 


and 

k.. = a positive integer, i =0, 1, 2, ..., n^-1. 

The corresponding expression for the output of the quantizer can 
be written as 

_ i o -(n-,-1 ) 

Y(z) = h(k b + k,b,z + k 0 b 0 z c + ... + k ,b ,z 
oo 11 2 2 n-|-l n-|-l 
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-n, -1 


(2n. -1 ) 


"1 "-I ' * / 

- kb z - k,b,z - ... - k -.b ,z + ...) 

oo 11 n-j - 1 n-j-1 


-n, -2n, ■, -(n,-l) 

h(l - z + z - •••)(k b 0 + k,b,z + ... + k„ ,b _,z ) 


'1 ' "1 


-n -1 n l -1 

= h(l + z ] ) { l k.b.z -1 ) 

i=0 1 1 

i 

n,-l 

n -j 1 • 

h'Zf ( I k.b.z - ) 
i=0 1 1 




where 

i 

h = level of quantization 

b^. = +1 or -1, depending on the phase of the input; i =0, 1, 



The variables b^ , i =0, 1, 2, ..., n^-1, determine the sign of 
the corresponding input and output pulses of the quantizer. The value 
of b.. is determined by the variable i and the phase shift <{> associated 
with the input cosine signal. 

The input to the sampler is written as 


e(t) = E cos (wt + <j)) 


(5-49) 
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and the input to the quantizer is 


e*(t) = E l cos (kwT + 4>}<5(t - kT) 
k=0 


For an oscillation of order n, we have 


uT = 


2ir 


or 


wT 


jl 

n. 


when n is even 


Thus, Eq. (5-50) becomes 


e*(t) = E l cos{^+ o,)6(t - kT) 
k=0 n l 


Let 


a. 

1 n l 


i =0, 1, 2, ...» n-j -1 


Then the variable b^ is given by 


b.j =■ SIGN[cos (a- + 4>)] 


With a given input amplitude E, an output of k^h at time a. 
requires the following relation to exist: 


(2k. - (2k. +1)h 

2 < ! E cos (a- + <f>) | < — — 2 


(5-50) 


(5-51) 


(5-52) 


(5-53) 


(5-54) 


(5-55) 


(5-56) 



Equation (5-56) implies that for a given even n, all combinations 
of k_.'s in the mode expression in Eq. (5-46) may not exist. Further, 
the range of <j> over which a particular mode exists will vary from mode 
to mode. These facts provide useful insight into the operation of a 
quantizer and play an important role in the implementation of any 

i 

scheme which is used to calculate the critical regions. 

In general, the range of <p which must be considered in describing 

function derivations is 2tt radians, but due to symmetry the required 

i 

variation of <f> can be reduced substantially. In the case of a quantizer 
and with n being even, the output waveform repeats after every 2n/n 
radians change in <J>. This fact was illustrated for n = 4 by the wave- 
forms of Fig. 5-4. Thus, it is sufficient to consider the range 
-ir/n <_<{)'_< n/n. 

The expression for the output of the quantizer, Eq. (5-48), is 
valid over that;./ range of 4> where the k^'s and b^'s do not change from 
their selected values. The expression in Eq. (5-56) provides the 
restriction on <f> because of possible variation in k. . To determine 
the influence of $ on b., it is necessary to determine minimum domains 
of constant b.. 1 s when <j> is varied. With reference to Fig. 5-7, the 

time of pulse occurrence of the quantizer input/output is according to 

* 

the sequence 


n 2 tt 4jt i 2tt (n-1 )2tt 

* n ’ n * ‘ ’ ‘ * n 5 * ’ ' ’ n 

or 


i 2 tt 


i = 0, 1 , 2, , 


n-1 
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With <j> = 0 as reference, the range of $ for constant is defined 
by the angle when any of the b^'s first change sign. Thus, the pulses 
located nearest to the tt/2 and 3tt/2 points must be considered; i.e., 
we need to determine those values of i for which |u/2 - 2m'/ n | and 
| 3 tt/ 2 - 2m/n| achieve a minimum. 

First, note that 


min 

i 


3tt 

2m‘ 

= min 

it . 2m o . 2iTn 

2 

n 

i 

2 n n 


min 

i 


ir , 2ir , 

- 2 + T (n 


- i) 




2frj 

n 


(5-57) 


i = 0, 1 , 2, . . . , n, j = 0, 1 , 2, . . . , n. 


Thus, it is sufficient to consider the minimum of J tt/2 - 2m /n| only. 
Two separate cases have to be considered: 


Ccue. a: nj e.ve.n 

When n^ is even, let 

n-j 3 2n r 4 (5-58) 


then 


n = 4n r (5-59) 

If we choose i = n^, = n/4, then 
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IT 

2-iTi 

IT 

7T 

2 ' 

n 

2 " 

2 


which is then the desired minimum. In this case, there is a pulse at 
or = tt/2 and also at cr = 3ir/2. The pulses nearest to these lie 
2-rr/n radians on either side; thus, the tr's will remain constant for 
the ranges 0 £ <j> £ 2-rr/n and -2-rr/n £ <p £ 0. Since the range of ' <J> 
variation is -tr/n <_<(»<_ rr/n, two ranges of <J> have to be considered 


Range 1 : * 0 £ <p £ 7T/n 

Range 2: -Tr/n £ 4> £ 0 (5-60) 

Ccue. b: Uj odd 

In this case the pulse nearest to tt/2 is given by i = (n^/2 ± 1/2), 

thus 


. 7T 

min j 
i 


2-rri 

m',' 


i . + 1, 

2 n [ 2 - 2‘ 


_ ir 2jT,n U 
2 " nV? 1 



7T 

n 


(5-61) 


Since the b.'s will remain unchanged for 0 £ <jj £ -rr/n and -rr/n £ (j) £ 0, 
the two ranges of <p are the same as in case a, Eq. (5-60). 

The ranges of <f> defined by Eq. (5-60) must therefore be considered 
separately when the expression for the quantizer output in Eq.(5-48) is 
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utilized. Although, the choice of the k/s further restricts the range 
of <}), this restriction will not be considered in advance, and will 
appear naturally when the critical regions are calculated. In fact, 
for some combination of the k^'s the range of <j> will be nonexistent. 

In order to derive a general expression for - 1 / N ( z ) when n is 
even, we write the expression for the input from Eq. (5-3) and the 
expression for the output from Eq. (5-48): 



(5-62) 




(5-64) 


(5-65) 
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Nl W 


Ezj sin u)T(cos (j) + j sin <j>) (z + 1 ) 


V 1 


(z 2 - 2z cos coT + l)h l k.b.z 1 

i=0 1 1 


n,-i 


Define 


C n 


1 

zj sin o)T(z + 1 ) 


n,-l „ . 

2 n l _1 

{z c - 2z cos coT + 1) l k.b.z 

i i=0 1 


I * i 

The term (z V 1) can be expanded as 


n, V 1 j(2k+l)7r/n 

z + 1 = n (z - e 
k=0 


1 


n,-2 


j it/ n-. -j-rr/n, 1 

= (z - e 1 )(z - e 1 ) n (z - e 

k=l 


j (2k+l )ir/n 


P V 2 j(2k+l)ir/n 1 

= (z - 2zcosu)T + 1) n (z - e 

k=l 


Substituting Eq. (5-68) into Eq. (5-67) yields 

i 


V 2 j (2k+l )w/n 


zj sin wT n (z - e 
k=l 


V 1 

j k.b.z 1 
i=0 1 1 


or 


(5-66) 


(5-67) 


(5-68) 


(5-69) 
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C 


n 


V 2 j{2k+l)Ti/n 1 

j sin ojT n (z - e ) 

k=l 


let 


'nR 


+ J'C 


nl 


where C nR and d n j are the real and imaginary parts of the constant C n 
The expression in Eq. (5-66) can now be written as 

- NriT = K (C nR + JC nI )(cos,(,+Js1n,|,) 


The constraints on the magnitude of E are 

(2k - 1 )h (2k . + 1 )h 

2 + 'l I E cos (<f> + ol . ) | < 2 

i = 0, 1 , 2, . . . , n^-1 ; k i f 0. 

Since E is a positive quantity, the expression in Eq. (5-73) can 
be written as 

(2k. - l)h 4 (2k + l)h 

2 < E l cos ( 4 > + a.) | < 2 

Let the range of <p be defined as 

4>1 < <p <_ <p 2 


where and ^ can take the values as specified in Eq. (5-60). 


(5-70) 


(5-71) 


(5-72) 


(5-73) 


(5-74) 


(5-75) 
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Define as 

„ *1 * *2 
*N ' 2 


The absolute value sign in Eq. (5-74) can be removed by use of 
new variable as 


a i = a i it < (^N + a i ) < 2tt or 0 < (^n + a i ) < f 

i 

and 


a;. = a.; 


TT 


if 


2 < ( a -j + 4>^) < 2 


3tt 


Equation (5-74) now becomes 

(2k - l)h (2k + l)h 

2 _< E cos (4> + a.j ) £ 2 

i'.' 

Thus, the maximum and minimum values of E can be written as 
(2k. + l)h 

r — L, 

max i 2 cos (4) + ) 

(2k - l)h 

min i 2 cos (<j> + aT ) 
i = 0, 1 , 2 n^-1 . 

The corresponding maximum and minimum values of -1/N(z) are 


1 

N 


(2k i + l)(C nR + jC nI )(cos 4) + j sin 4>) 
— 2 cos (<j> + ) 


(5-76) 


(5-77) 


(5-78) 


(5-79) 

(5-80) 


max i 


(5-81) 
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1 

N 


min i 


(2k . - 1 )(C nR + jC nI )(cos <j> + j sin <*>) 
2 cos (* + ai ) 


In order to eliminate * from Eqs. (5-81) and (5-82) consider 


Eq. (5-81) first; expanding the denominator and numerator yields 


1 

N 


max i 


(2k. + 1 )[(C nR cos <p - C nI sin*) + j(C nR sin<() + C nI 
2(cos * cos aj - sin * sin ) 


Let > 


1 

N 

max i 


I 

= x + jy 


By comparing Eqs. (5-83) and (5-84), we have 


( 2k i + C nR cos * " C nl sin<<> 

2 cos <p cos a! - sin * sin aT 

K + C nR “ C nl tan(|) 
l cos a! - sin ai tan * 


where 


< = (2k. + 1 )/2 

i 


and 


y = 


C nR s1n»^ C nl cos» 
i cos * cos al - sin * sin a.j 


, C D tan (p + C T 
= K + nR J nl 
i cos aT - sin aT tan * 

Equation (5-85) is solved for tan* to yield 


(5-82) 


COS <j>) 


(5-83) 


(5-84) 


(5-85) 


(5-86) 


(5-87) 



108 


tan = 


K + C D 
l nR 


K + C T 
l nl 


x cos a I 
x sin a! 


(5-88) 


If Eq. (5-88) is substituted into Eq. (5-87), we have 


Y = 


C nR 

- !/ + 

+ ) 
K. C n - x cos a! 

^ nR l 

+ c T 

nl 

K + C T - x sin a'. 
i nl i 

" "i 

cos - sinaj 

i 

C nR ‘ x cos “1 

K + C T - x s i n a '. 
i nl i 


(5-89) 


After some algebraic manipulation, Eq. (5-89) yields 


2 2 

r c + r 

u nR nl 


'i C nI cosa| - Casino: 


C ~ cos a I + C T sin a 1 . 
nR i nl i 

C T cos a! - C D sin a'. 
nl i nR ^ 


(5-90) 


Equation (5-90) which is obtained from Eq. (5-81) by eliminating 
<f) shows that the.' -1/N|max i locus in the complex plane is a straight 
line. Similarly, when <J> is eliminated from the -1 /N jmi n i expression, 
Eq . (5-82) , we have 


K i 


°nR 

C „I cos “1 


+ C 


nl 

C nR sin “i 


C nR cos “i + C nl sln a i 
C nI c0Sa i - C nR sina i * 


4 

which is also a straight line in the complex plane, and 


(2k. - 1) 


(5-91) 


(5-92) 


Note that the lines defined by Eqs. (5-90) and (5-91) are parallel 
but have different intercepts on the imaginary axis. Thus, each magnitude 
constraint of Eqs. (5-79) and (5-80) on the input expression of 



Eq. (5-62) generates two parallel lines which represent the maximum 
and minimum values of -1/N(z), respectively. 

Since there are n^ magnitude constraints, Eq. (5-73), each 
output expression generates 2n^ lines, where n-| lines are the maximum 
values of -1/N(z), and n^ lines are the minimum values of -1/N(z). 

When n is even, the number of <j> ranges to be considered is two, 
as in Eq. (5-60), then Eq. (5-63) gives two output expressions which 
have to be considered. The difference between these two output 

t 

expressions would be that some of the b . 1 s are different due to the 

i * 

different <j> rahge. Therefore, each mode as defined by Eq. (5-46) 
yields the following number of lines: 

n^ = 4n-| (5 

where n^ represents the number of lines generated in the complex 
plane. 

Among the lines in Eq. (5-93), 2n^ represent maximum conditions 
for -1/N(z), and 2n-| represent minimum conditions. Also, each maximum 
line is parallel to the corresponding minimum line. 

The critical region for -1/N(z) is thus the region enclosed by the 
n. lines and satisfying the maximum and minimum requirements associated 

L 4 

with each line. In some cases if the choice of the k_.‘s in Eq. (5-46) 
is such that there may be no region enclosed, then that particular mode 
will not exist. 

The procedure for determining the critical region for the quantizer 
describing function when n is even, can be summarized as follows: 

1. Select an even integral n. 



no 


2. Select a mode by choosing the k^'s in Eq. (5-46). 

3. There are two ranges of <j>, 0 £ <j> £ u/n and -Tr/n _< 4> <_ 0. 

For each range of <f> determine the b^'s and thus the output, 

Eq. (5-63). 

4. Calculate the coefficient C n as in Eq. (5-70) for each range 
of 4>. 

5. Determine the a.'s, Eq. (5-54), and aj's, Eq. (5-77), for the 
chosen n. 

. i 

6. For each and each range of <f> , calculate the slope and 
intercepts of the n^ lines, Eqs. (5-90) and (5-91) which define 
the maximum and minimum. values of -1/N(z). 

7. Determine the region enclosed by these lines. 

This is the critical region. 

As n is increased, the number of equations defined by Eq. (5-93) 
quickly rises a^pd it is useful to consider utilizing a digital computer 
for determining the critical regions. The following table illustrates the 
number of lines associated with various even values of n > 4. 


n 

ranges of <j> 
(deqrees) 

no. of lines 
per <J> range 
(2n-| ) 

total number 
of lines 

( 4n l ) 


i 



4 

0 £ 4> £ 45 

8 

16 


-45 < <j> £ 0 



6 

0 £ <j> < 30 

12 

24 


-30 <_ <j> <_ 0 



8 

0 < (j) < 22 . 5 

16 

32 


-22.5 £ <}> < 0 



10 

0 £ 4> £ 18 

20 

40 


-18 £ <p £ 0 



20 

0 £ <f> £ 9 

40 

80 
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B. Odd n 

When n is odd, the input pulse train to the quantizer can assume 
a maximum of n different values. The mode of oscillation is thus 
characterized by 


A (k Q » k-j , k^» •••> k n _i ) 


(5-94) 


The output of the quantizer is written as 


Y(z) = h(k b + k,b,z~^ + k 0 b 0 z"^ + 
Voo 11 2 2 


+ k ,b 

n-1 n-1 


+ k b z~ n + k,b 1 z - ^ n+ ^ + ...) 
oo 11 


= h(l + z" n + z _2n + ...)(k 0 b 0 + k^z' 1 + ... + k_ 1 b_ 1 z" (n “ 1) ) • 


n-1 n-T 


h( V“ z' n )(V k.b z’ 1 ) 
1=0 


_ n-1 

hz n ( l k.b z" 1 ) 
i=0 1 1 

z n - 1 


n-1 

h Y< k.b.z 
it 


n-i 


i=0 


z n - 1 


(5-95) 


where h and b.. have the same meaning as in the even n case. 
The input to the sampler is again written as 


e(t) = E cos (wt + <J>) 


(5-96) 
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and the input to the quantizer is 

co 

e*(t) = E l cos (kioT + <j>)6(t - kT) 
k=0 


= E l cos (^p- + 4>)6(t - kT) (5-97) 

k=0 n 

(5-98) 
(5-99) 

T 


(5-100) 

As in the case of even n, Eq. (5-100) restricts the possible 
combinations of k^'s in the mode expression of Eq. (5-94). 

In order to^determine the range of <J>, consider the waveforms 
in Fig. 5-8. In this figure n = 3 and the waveforms repeat every 120 
degrees. In general, for odd n, the output waveform will repeat 
when <p is changed by 2Tr/n radians. Thus, the range of <j> is chosen as 
-u/n <_ <j> £ Tr/n . 

It is now necessary to divide the chosen range of <j> into smaller 
ranges depending on the variation of b. 's in the chosen full range. 


Let 


a i 


2ni 


i =0, 1 , 2, . . . , n-1 


then, b. is given by 


b.. = SIGN[cos (a^ + 4))] 


The input and output of the quantizer are constrained by the 
following relation: 


(2k. - l)h (2k. + l)h 

2 < l Ecos ( a i + 4*)| < 2 
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The effect of k. 's on restricting the range of <j> will be accounted 
for when the limits on the input magnitude as described by Eq. (5-100) 
are considered. To determine the influence of b^'s on <j>, we again 
calculate the domains of constant b^ when $ is varied in the chosen 
range. The pulse occurrence timing is written as 

n (n-1 )2 tt 

I 

i = 0, 1 , 2, . . . , n-1 . 

Taking <J> = 0 as reference, it is necessary to determine the values 
of i for which |tt/2 - 27ri/n| becomes a minimum. The result of Eq. 
(5-57) is still valid and allows us to ignore the 3tt/2 point while 
considering the jt/2 point only. 

The minimum of |tt/2 - 27ri/n| will be achieved by one of the 
two pulse occurrences on either side of the ir/2 points (when n is odd, 
no pulse occurs at tt/2 exactly). Thus, the minimum will be given by 



or 

2iri 
n * 


t 

which ever sign makes i an integer. Thus, when n = 3, the + sign is 
used and when n = 5 the - sign is used, i = 1 in each case. Thus, 


IT 

27Ti 


TT 

2tt 

[n + ll 

2 ' 

n 


2 " 

n , 

[4 - 4j 



_7T_ 

2n 


(5-101) 
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The chosen range of <j>, -ir/n £ 4> £ ir/n has to be subdivided into 
sections given by Eq. (5-101), over which the b^'s are constant. 

Four ranges of <j> must therefore be considered. 


Range 1 : ir/2n £ <{> £ ir/n 

Range 2: 0 <_ <j> <_ n/2n 

Range 3: -Tr/2n £ <f) £ 0 

Range 4: -ir/n £ <j> £ -7r/2n (5-102) 

t 

A general expression for -1/N(z) with odd n is now derived. The in- 
put signal will change when <J> is varied beyond its range as defined by Eq. 
(5-102). Accordingly, each range in Eq. (5-102) will have to be 
calculated separately. 

The input and the output of the quantizer are written as 


E(z) 


Ezj sin uiT(cos <t> + j sin <J>) 
2 

z - 2z cos u)T + 1 


(5-103) 


h V k.b.z 0 " 1 ' 

Y(z) = — 1 Z.9 — (5-104) 

z n + 1 

respectively. 

The negative inverse of the describing function is 

1 = . Hz 

N(z) Y[z 



Ezj sin ojT(cqs 4) + j sin <t>)(z n + 1 ) 
? n-1 

(z^ - 2zcosa)T + l)h l k.b.z 11 " 1 

i=0 1 1 


(5-105) 


Define 
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C n = - 


zj sin o)T(z + 1 ) 


? n-l • 

{z c - 2z cos ojT + 1)( l k.b.z n " } 

1=0 1 1 


In Eq. (5-106), the term (z 11 - 1) is written as 


n-l 


(z n - 1) = n (z - e^ 2TTl/n ) 
i=0 


n 2 

- (z - l)(z 2 - 2zcoscoT + 1) IT (z - e j27ri/,n ) 

i=2 


Substituting Eq. (5-107) into Eq. (5-106) yields 


C n 


n 2 

zj sin wT(z - 1) n (z - 
i=2 


n-l 

l V 

i=0 1 


n-i 


n 2 

j s i n wT ( z - T) n (z - e j2 ' rT ' i / n ) 

- i^2 

l k i-l b i-l zn_i 


Dividing Eq.,. (5-108) into real and imaginary components, we 

4 


wri te 


C n C nR + jC nI 


The expression in Eq. (5-105) then becomes 


NTD- = K (C nR + jC nI> (cos< l , + js1n, '' ) 


(5-106) 


(5-107) 


(5-108) 


(5-109) 


(5-110) 
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The magnitude of E is constrained by 

(2k. - 1 )h (2k + l)h 

2 1 |E cos (<() + cu ) | < ^ 

i = 0, 1 , 2, . . . , n-1 , k. j 4 ' 0. 

Let the range of <}> under consideration be 

1 <P < <t>2 

i 

where and tp? can take on values as specified in Eq. (5-102), we 
define •• 

^1 + 

*N = — 2— 

Select a! such that 

a i = a i l f T" < + a i^ < 2it 

or 0 < ((j> N + a.j ) < tt/2 

a \ = - it if < (<^ + ou) < 3 tt/2 

Then, Eq. (5-111) can be written as 

(2k. - l)h (2k, + l)h 

2 £ E cos (<}> + al ) _< — ■ — 2 

The maximum and minimum values of E from Eq. (5-115) become 


(5-111) 


(5-112) 


(5-113) 


(5-114) 


(5-115) 
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( 2 k +l)h 

r _ I 

max i " 2 cos ( <j> + a] ) 

( 2 k. - l)h 

c = !_ 

min i 2 cos (<j> + a- ) 

The corresponding maximum and minimum values of -1/N(z) are 

1 (2k. + 1 )(C nR + jC pI )(cos ()> + j sin 4 ,) 

N max( ' 2 “ s U + “t ) 

1 (2k i - 1 )(C nR + jC nI )(cos 4 > + j sin <j>) 

’ W min i = 2 cos + “i > 

Equations (5-118) and (5-119) are similar to Eqs. (5-81) and (5-82), 
respectively, obtained earlier for even n. Thus, as in that case, it 
is possible to eliminate 4 ) and determine the straight lines defined 
by the locus ofi'.each expression which 4 > is varied. The result is 
identical to Eqs. (5-90) through (5-92). 

Each magnitude constrai nt in Eq. (5-115) again generates two 
parallel straight lines, one representing the maximum value of - 1 / N ( z ) 
and the other representing the minimum value of -1/N(z). 

With n magnitude constraints and four ranges of 4 >, the total 

4 

number of straight lines in the complex plane is 

n L = 8 n (5-120) 

As before, the critical region for -1/N(z) is the region enclosed 
by the n^ lines and simultaneously satisfying the maximum and minimum 
requirements associated with each line. Also, if the choice of k^'s 
is such that no region may be enclosed, that particular mode will not exist. 


(5-116) 

(5-117) 


(5-118) 

(5-119) 
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The procedure for determining the critical region for the quantizer 
describing function when n is odd is summarized below: 

1. Select an odd integral n. 

2. Select a mode by choosing the k^'s in Eq. (5-94). 

3. Four ranges of <f> can exist, as defined by Eq. (5-102). 

Determine the b^'s for each range. 

4. Calculate coefficients C n as given by Eq. (5-108) for each 
range of <j>. 

15. Determine the oit's, Eq. (5-114). 

6. For each at and each range of <j> calculate the two lines of 
Eqs. (5-118) and (5-119) which define the maximum and minimum 
values of -1/N(z). 

7. Determine the critical region by finding the region enclosed 
by all the 1 ines. 

The following table indicates the number of lines associated with 
several values of odd n: 


n 


Ranges of <j> 
(degrees) 


no. of lines 
per <j> range 
<2n) 


total number 
of 1 i nes 
(8n) 


3 0 j< _< 30 , 30 £ 4) £ 60 6 

-60 £ <|) £ -30 , -30 £ <J> £ 0 

5 0 £ (p £ 18 , 18 _< 4> _< 36 10 

-36 £ <J> £ -18 , -18 £ <j> £ 0 

7 0 < <p £ 12.9 , 12.9 < <p < 25.8 14 

-25.8 < <p < -12.9, -12.9 < <fi < 0 


48 

80 
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5-5. Computer Implementation and Results with the General Quantizer 
Describing Function 

The DDF derived in the previous sections has been programmed on a 
digital computer. Input data consists of the order of oscillation n 
and the modes which are to be investigated. The program calculates 
the cj) ranges, the aVs, the b^'s and generates all the possible lines 
for the given mode. These lines are then plotted by hand and the 
critical region determined. 

As examples of the utility of the program, two values of n, 
n = 4 (even case) and n = 3 (odd case) are selected. 

A. Results for n = 4 

Figure 5-9 shows a typical printed output from the program for 
one mode, A = (1, 2). In Fig. 5-9 the following information is 
provided: 

1. N = the order of oscillation n 

N1 = n/2 for even n and N1 = n for odd n 

NPR = no. of <p ranges 
NLS = NPRxNI 

NL ' = no. of 1 ines n. 

L n (n - 2) 

N1NT = no. of intersections of the lines, 2 

2. The cp range is printed nextasPHI RANGE 

3. The ct^.'s are written as ALPHA VALUES 

4. Sequence of b Q k o , b-jk-j, ... is written as SEQUENCE OF 

KO, K1 , etc. 

5. C n £ and C n j are printed next 
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DISCRETE DESCRIBING FUNCTION FDR. A QUANTIZER NONLINEARITY 

N= 4 Nl= £ NPF- £ . NLS= 4 Nl = S NINT= £4 

PHI RANGE IS-45.00 TO 0.00 

ALPHA VALUES ARE 0.000E-01 9.000E+01 

SEQUENCE OF K 0 f K 1 »...ETC IS 1. £. 

• C N - R E A L = - £ . 0 0 0 0 0 E - 0 1 C N - 1 M A G = - 4 . 0 0 0 0 0 E - 0 1 
THE LINES ARE GIVEN BY Y=X»SLDPE+CRDSS 

SLOPE CROSS 1 CROSSE C 

-5.0 0 0 0 0E- 01 -7 . 49999E- 0 1 -£ . 5 0 0 0 0E- 0 1 

2 . 0 0 0 0 1 E + 0 0 £.50 0 0 1 E + 0 0 1 . 5 0 0 0 0 E + 0 0 

PHI RANGE IS O.UO TO 45.00 
ALPHA VALUES ARE 0 . 0 0 0E- 0 1 -9 . 0 0 0c+ 0 1 . 

SEQUENCE DF KO jKI > . . .ETC IS 1. -2. 

CN-RE AL=-£ . 0 0 0 0 0E- 0 1 CN- I MAG= 4 . 0 0 0 0 0E- 0 1 
THE LINES ARE GIVEN BY Y=K*SLOPE+CROSS 

SLOPE CROSS 1 CROSSE 

5 . 0 0 0 0 0 E - 0 1 7.5 0 0 0 1 E - 0 1 £ . 5 0 0 0 0 E - 0 1 

- £ . 0 0 0 0 1 E + 0 0 - £ . 5 0 0 0 1 E + 0 0 - 1 . 5 0 0 0 0 E + 0 0 


Figure 5-9. 
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The lines for this 4* range are then written, where SLOPE 
refers to the slope of the maximum and minimum lines. 

CROSS! is the intercept of the maximum line and CR0SS2 is 
the same for the minimum line. The lines are written in the 
same order as the a!'s, or the occurrence of the pulses. 

7. All information is repeated for the other <j> range. 

The lines are plotted and the critical region is determined. 

Figure 5-10 shows the plots for several different modes with n = 4. 

In this case the critical regions are symmetrically located around the 
(-1, 0) point in the complex plane. Also, as k Q , k-j increase the 
critical region tends to shrink towards the (-1, 0) point. 

B. Results for n = 3 

With n = 3 the mode of oscillation is defined by three different 
variables, k Q , k-j and Also, in this case, many of the combinations 
of k^ 's cannot exist. Figure 5-11 shows the critical regions obtained 
with the computer program for A = (1, 0, 0) and A = (1, 1, 1). The 
critical region was found not to exist for the modes A = (2, 0, 0), 

A = (1, 2, 2), A = (1, 2, 3). 

Results for other modes and other values of n can be obtained 
similarly. However, as n gets larger, the number of lines increases 
substantially and the procedure for determining the critical region 
gets more involved. It would be desirable if the computer program were 
able to calculate these regions itself from the straight lines. In 
order to expand the scope of the program to include this calculation, 
several special cases, which can occur have to be resolved first: 









1. In certain modes, the duplication of the lines results in 
many parallel lines, and their intersection calculations 

can cause singularities. Such duplicate lines and repetitions 
should be recognized and eliminated, thus removing any 
singularities as well as reducing the computation. 

2. In certain modes, the straight line can be a line with infinite 
slope. This happens when certain combinations of C^, C j and 
a. 's yield zero denominator in the slope expression, Eq. (5-90). 
This condition should be recognized and the line should be 
redefined as x = f(c, a), with no occurrence of y in the line 
equation. 

3. Sometimes the straight line has zero slope and this must 
also be recognized and eliminated. The line must be re- 
defined using y only. 

These comments can provide a basis for extending the results 
presented here and using them in a more general situation. 
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